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Manipulating a single adsorbed DNA for a critical endpoint 
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Abstract. - We show the existence of a critical endpoint in the phase diagram of unzipping of 
an adsorbed double-stranded (ds) polymer like DNA. The competition of base pairing, adsorption 
and stretching by an external force leads to the critical end point. From exact results, the location 
of the critical end point is determined and its classical nature established. 
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A critical end point (CEP) occurs when a continu- 
ous transition line terminates on a first order transition 
ine. It's specialty: even though it is surrounded by three 
phases, still there is no three phase coexistence, but, in- 
stead, a scale-free critical phase coexists with a noncriti- 
cal phase [1,2]. For comparison, a triple point, also sur- 
rounded by three phases, would show three phase coex- 
istence. A CEP is expected to occur in various mixtures 
and ferroelectrics and in vortex lattice [1—3]. We show 
here that a different way of locating and studying a CEP 
is through single molecular manipulations of an adsorbed 
DNA. 

The melting of DNA is known to be a crucial step 
in many biological processes [4]. The double stranded 
DNA(dsDNA) is a bound state of two polymers or strands 
held together by hydrogen bonds of base pairs. The phe- 
nomenon of cooperative breaking of the base-pairings ther- 
mally or otherwise is melting. The recognition of force 
as a thermodynamic variable for this process has helped 
in completing the phase transition picture of dsDNA [5] . 
Even though the nature of the thermal denaturation of 
DNA remains a puzzle, the force induced unzipping tran- 
sition is theoretically well-settled [5-18]. In addition, the 
force induced unzipping transition has emerged as a pos- 
sible scenario for opening of DNA [19] with the replication 
Y-fork as the junction of two phase coexistence, the zipped 
and the unzipped DNA. A thermodynamic description of 
DNA would entail two conjugate ensembles of fixed force 
and fixed distance. These two ensembles are important in 
both theoretical and experimental situations [8,20-22]. 
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The Melting and unzipping of DNA are generally con- 
sidered in the free environment of bulk solutions, but often 
the presence of interacting surfaces cannot be ignored. In 
vivo, during replication, DNA gets attached to the mem- 
brane but otherwise it remains away ( "desorbed" ) from 
the membrane. The protein-induced membrane-DNA at- 
tachment is used in the replication process and cell divi- 
sion [23]. In gene therapy, targeted delivery is achieved 
by taking advantage of adsorption-desorption of DNA on 
cationic liposomes [24,25]. That metallic (e.g., gold), semi- 
conducting (e.g., silicon) or insulating (e.g., mica) surfaces 
can also adsorb DNA has opened up the possibility of 
biosensors for fast and precise detection of DNA in samples 
like hair, blood etc. In all these cases, the surface-DNA 
interaction depends (and hence tunable) on the nature of 
the surface, fluctuation of the surface as for fluid mem- 
branes, ionic concentration of the environment, nature of 
hydrophobicity and van der Waal interactions. A well- 
studied system in recent years is DNA on gold where the 
DNA can be attached to the surface with a thiol group and 
a small linker [26-28] . The model setup we are considering 
is similar to the case of DNA on a gold plate and is shown 
in Fig. 1. A force pulls only one of the two strands of a 
DNA which can adsorb on an attractive surface. The DNA 
as usual has the base pairing energy. We here treat the 
surface-DNA interaction as an additional parameter in the 
problem. It transpires that for certain ranges of the attrac- 
tion with surface there would be a competition between 
adsorption, unzipping, and melting. These three processes 
can lead to a CEP. The essential feature that holds the key 
for the CEP is that the force-induced unzipping is strictly 
first order but the adsorption-desorption transition of a 



Rajeev Kapri and Somcndra M. Bhattacharjec 



polymer from a surface is continuous [29-31]. 

Previous studies of DNA unzipping showed that the lat- 
tice model preserves, even in two dimensions, the basic 
results of DNA unzipping including the first order nature 
of the phase transition and the existence of a reentrant 
region [7]. The generic arguments of these studies also 
showed that the choice of the lattice is not crucial. For 
the problem at hand, we have also done Monte Carlo sim- 
ulations in 2 + 1 dimensions (on a cubic lattice) and find 
that the force-distance isotherms are qualitatively similar 
to the isotherms obtained in 1 + 1 dimensions (see be- 
low). Consequently we focus on the results obtained from 
analysis of exact results in 1+1 dimensions. 

We model the DNA by two directed self avoiding walks 
on a D = 1 + 1 dimensional square lattice. See Fig. 1. The 
walks, labeled 1 and 2, starting from the origin, are di- 
rected along the diagonal of the square (z direction). The 
walks are not allowed to cross each other but whenever 
they meet (i.e. x\{z) = £2(2)) there is a gain in energy 
—e b (e b > 0) for every contact. This is the base pairing. At 
the diagonal (x = 0) there is an impenetrable attractive 
surface, with an energy — e w (e w > 0), which favours the 
adsorption of the DNA. In 1 + 1 dimensions, the surface 
is a line passing through the diagonal of a square lattice, 
and only one of the strands can get adsorbed on it (i.e. 
xi =0), since the two strands of the DNA cannot cross 
each other. One end of the DNA is always kept anchored 
at the origin. We apply an external force g, along the 
transverse direction (x-dircction) on the free end of one 
of the strands of the DNA. The other strand is left free. 
Henceforth, the strand which is left free is called the "free 
strand" and the strand, on which the external force acts is 
called the "pulled strand" . The endpoint positions Xi (N) 
is to be shortened to Xj. In 2 + 1 dimensions, the surface is 
a plane passing through the diagonal of a cubic lattice [31]. 
Unlike the 1 + 1 dimensional case, both the strands of the 
DNA can get adsorbed on the surface and still satisfy the 
non-crossing constraint on the plane (y-direction). 

The two energies independently give us two special tem- 
peratures: (i) T w the temperature for desorption of the 
DNA from the surface, and (ii) T m the melting tempera- 
ture of dsDNA. In absence of a surface, the melting tem- 
perature is given by k B T m = e 6 /ln(4/3) [7-9]. For ssDNA, 
fcfiT w = e w / ln2 [29]. Thermal fluctuations create bubbles 
in the dsDNA changing its effective elastic behaviour with 
concomitant rise in the desorption temperature. We con- 
sider energies such that T w < T m and for numerical results 
we choose units so that ks = 1, and 66 = 1. This conven- 
tion is also adopted in the following discussion, unless we 
want to show a general formula or show the dependence 
on e b . 

There are four distinct phases differentiated by (xj) for 
% = 1,2 as N — > 00, ((...) denotes thermal averaging), (i) 
Z a : zipped DNA adsorbed on the surface with (xi)/N — ► 
for i = 1,2. (ii) Z d : zipped DNA desorbed from the sur- 
face with (xi)/N = 0(1), for i = 1,2. (Hi) U a d: un- 
zipped DNA with the free strand adsorbed on the surface 



and the pulled strand is stretched in the direction of the 
force. This phase is characterized by (xi)/N — ► but 
(x 2 )/N = 0(1). (iv) XJ dd : unzipped DNA with both the 
strands desorbed from the surface with (x2)/N = O(l) 
and (xi) /\/~N — O(l). The adsorbed or the zipped phases 
may also be characterized by the fraction of monomers 
in contact, e.g., $ w the fraction of polymers in contact 
with the wall and <J>b the fraction of bound base pairs. 
In a zipped phase $b / while for the adsorbed phase 

One can do a zero temperature (T = 0) analysis of 
the problem, keeping e b constant. The energies of the 
three phases, namely Z a , Zd, and U a( j are respectively 
given by E Za = -N(e w /2 + e 6 ), E Zd = -N(g + e 6 ), and 
£x, ad = -N(e w /2 + g). For e b < e w < 2e b , the phase 
U a d is always unstable (i.e. it has higher energy). The 
transition from phase Z a to phase Zd occurs at g = e w /2. 
On the other hand, for e w = 2e b , there is a degeneracy for 
Zd and U a d which occurs at g = e w /2. But, for e w > 2e b , 
Zd phase is no longer favourable. The change in stability 
of phases Zd and U a d CIS Cyj IS tuned is an indication that 
these phases could be stabilized by entropy at intermediate 
temperatures under appropriate conditions. For e w = 0, 
we can only have Zd and Udd, while, for e w > 0, as the 
T = analysis shows, phase Z a must exist. For very high 
temperatures T T w , T m , the stable phase is Udd- When 
e w = 00, the free strand, which remains adsorbed on the 
surface at all temperatures, acts like a zig-zag hard-wall. 
In this case, only phases Z a and U a d survive in the phase 
diagram. Consequently, there is a continuous evolution 
of the phase diagram for DNA £LS Cyj IS changed. These 
phases can therefore be represented in a 3-dimensional 
g-T-e w phase diagram. We show the cross-sections (g-T 
plane) of this phase diagram for various e w in Fig. 2. 

We like to point out the generic nature of the results, 
especially the existence of the phases and the nature of 
the phase transitions. Similarly, the choice of a straight 
wall is not a restriction; a zig-zag wall with all monomers 
getting absorbed also show similar behaviour [33]. For an 
experimental realization, DNA tethered to a gold or mica 
surface looks promising. For example, mica-DNA interac- 
tion can be fine tuned by NaCl and MgCl2 over a broad 
range of 0.02ev per bp to 0.35ev per bp [28]. For gold 1 
kbp DNA adsorption has also been studied [26]. Unzip- 
ping force force measurements via atomic force microscopy 
for such finite chains with proper finite size analysis [12] 
could verify the theoretical results presented here. 

Let D n (x\ 1 X2) be the partition function (temperature 
dependence not shown explicitly) of a dsDNA in the fixed 
distance ensemble where nth monomers of the strands are 
at positions x\ (free strand) and xi (X2 > xi) (pulled 
strand) respectively from the wall. D Il (xi,X2) satisfies 



1 There is a fifth phase U aa in 2 + 1 dimensions where DNA melts 
but both strands remain adsorbed on the surface. See, e.g.. [32] 
for experimental determination of such melting temperatures of ad- 
sorbed oligomers on gold. 
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the recursion relation ( x 2 > x\ > 0) 

D n+ i(xi,X 2 ) = ^ D n(xi+i,x 2 +j) 
i,j=±l 

x [1 + W6 X1 , ][1+B6 X1 , X2 ], (la) 

where, 

W = {ef )tv -l) 1 B={e f3t » -1), and/3=l/fc s T. (lb) 

The initial condition Dq(xi, x 2 ) = S Xlt oS X2t Q. The canoni- 
cal partition function with an external force g at the end 
of the pulled strand is then obtained by summing over all 
the allowed configurations of the DNA of length TV on the 
lattice. 



Z N (j3,g)= ^2 D N {x\,x 2 ) e 



0gx 2 



(2) 



From the partition function we calculate the endpoint av- 
erages (xi) and (x 2 ). The appropriate response function 
is the isothermal extensibility, which can be expressed in 
terms of fluctuations of the position of the end monomer 



X 



d(x 2 ) 



x 2 > - {x 2 ) 



(3) 



The N (length) dependence and finite size scaling of these 
quantities would be utilized to identify the phases and the 
phase transitions with the help of the Bhattacharjce-Scno 
data collapse program [34]. (See Fig. 3(a) as an example of 
the data collapse of x(#, T).) The nature of the transition 
(first or second order) is inferred from the values of the 
relevant exponents [12]. The phase diagrams are obtained 
by the repeated use of finite size scaling. 

In Fig. 4(a), we have shown the force-distance isotherms 
for e w = 1.8 at two different temperatures T = 0.5 and 1.5 
for the chain of length N = 1000. In both cases, we start 
from the ground state at g = 0, an adsorbed DNA on the 
surface. The phases can be identified by the extensivity 
of (xi) and/or (X2). At T — 0.5, there is a critical force, 
g s , at which the DNA gets unzipped from the surface but 
remains double stranded. We call this as "transition Sz". 
But at T = 1.5, we see the sequence Z a ^==> Zd U a d, 
an additional transition (to be called "transition Uz") at 
g = g u . The isotherms are obtained at T = 1.5 by two dif- 
ferent methods with results comparing nicely. The lines 
are from the exact transfer matrix based on Eq. (la), 
whereas the bigger symbols (squares and circles) are ob- 
tained by performing Monte Carlo simulations for longer 
chains using the multiple histogram technique [35]. The 
details of Monte Carlo simulation will be discussed else- 
where. The estimates, so obtained, are shown by the upper 
and the lower triangles for the free and the pulled strand 
respectively in Fig. 4(a). An isotherm for a 2 + 1 dimen- 
sional case is also shown in the Fig. 4(b) and is similar 
to the 2-dimcnsional CM SO , MS already mentioned. 

Figure 2(a) shows the phase diagram for e w = 1 as a rep- 
resentative in < e w < eb = 1. It contains three phases, 



namely Z a , Zd and Udd- The phase boundary separating 
phase Z a from phase Zd is shown by circles. The mini- 
mum temperature, T u , above which the unzipping of the 
dsDNA to two single strands takes place is at T u = e&/ In 2 
same as for e w = case. For low values of g, the pulled 
strand is not necessarily straight. If we ignore the effect 
of the wall the phase boundary can be calculated exactly 
from the above recursion relation as 



9u{T) 



2 



In 



2e _/3eb 



1 - 2e-P*> 



(4) 



Below T u , the DNA remains double stranded for any value 
of force g. That the effect of the wall is negligible is borne 
out by the excellent agreement [33] of the phase boundary 
for e w < 1 (Fig. 2a). 

The situation of interest is 1 < e w < 2. In this case, we 
have all the four phases in the phase diagram, including 
the T = unstable phase U a d- With e w > e& = 1, pair 
breaking, as opposed to desorption, would play an im- 
portant role at low temperatures. As a result, the phase 
boundary separating phase Z a from phase Zd, loses its 
monotonicity seen in the e w = 1 case and a thin slice of 
reentrance starts appearing in the phase diagram at inter- 
mediate temperatures. Such an intermediate reentrance 
was also found in Ref. [8], in a different context of a ds- 
DNA with a force at an interior point. We have shown the 
phase diagram for e w = 1.8 in Fig. 2(b), as a representa- 
tive of this regime with the reentrance shown in the inset. 
Apart from this feature, there is a region in the phase 
diagram which involves three phases, namely Z a , Zd and 
U a d- The transitions from phase Zd to phases U a d, and 
Udd are of first order, whereas, the transition from phase 
U a d to phase Udd is second order. The critical line is at 
T c = e w /lxi2 for all g. This is the temperature at which 
the adsorbed free strand in phase U a d desorbs from the 
surface. If we ignore the effect of attraction of the wall on 
phase Zd, the first order phase boundary can be obtained 
by equating the free energies of the two phases obtained 
from Eq. (la). For the Zd to U a d transition we then get 



m{T) = ^ In 



5 /3e b _ f e -/3e w _ g-2/3£„) 1 /2 



(5) 



while the Zd to Udd transition is still given by Eq. (4). 
Of course Eq. (5) will not be appropriate in the temper- 
ature region where there is reentrance in the desorption 
boundary. This is seen when compared with the numeri- 
cal results (Fig. 2b). However the large force asymptote is 
given correctly, especially the dependence of T u on e w as 



k H T u In 



1 - y/1 -4exp(-e 6 /fc B T u ) 



(6) 



This expression has the correct limits and matches with 
the numerics as shown in Fig. 3(b), justifying a posteriori 
the neglect of the effect of the attractive wall in these 
conditions. 
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The free strand desorption critical line terminates on 
the first order boundary for phase Zd separating it from 
phases U a( j and Udd- This point of intersection is the 
CEP. Note that g e (T) and g u (T) meet at T c = e w /\n2 
with same slope as they should for a CEP. This point is 
shown in Fig. 2(b), by C c (filled square). For e w = 1.8, 
C e is at g e = 0.742618..., T e = 2.59685... (with e b = 1). 
The CEP appears in the phase diagram only for e w > 1, 
and shifts towards the melting point as e w is increased. 
That the two curves meet at CEP with the same slope is 
a confirmation of its nature. At a triple point the angle 
between any two phase boundaries is strictly less than 2w. 

The phase diagram for e w = 2 has a new feature that 
Zd just becomes unstable at T = 0. For e w > 2, a triple 
point appears in the phase diagram where Z a , Zd and 
U a d coexist. The CEP still persists but Zd is now sta- 
bilized by entropy. With the increase of e w beyond 2, 
the region representing phase Zd shrinks rapidly and both 
the triple point and the CEP shift towards higher temper- 
atures and disappear independently from the phase dia- 
gram. The crossing of the thermal desorption and thermal 
melting temperature introduce new complications. These 
will be discussed elsewhere. 

In Fig. 4(c), we have plotted the scaled distances 
( X i,2)/N, of the end monomers of both the strands from 
the surface as a function of temperature T at an external 
applied force g — 0.925. This particular value of the force 
lies in a small region which allows us to see all the pos- 
sible phases, including the reentrance between the phase 
Z a and the phase Zd (see Fig. 2(b)). Just by increasing 
the temperature, the DNA can be made to go through the 
sequence 

Zd Z a <S=>- Zd U a d Udd- 

In the last phase, the free strand of the DNA desorbs from 
the surface and stays at a distance of y/N from the surface 
(not visible in this scale). To make it visible, we have plot- 
ted, in the inset, the scaled separation, (xi)/vN, of the 
end monomer of free strand from the surface as a function 
of T for DNA of various lengths. The plot confirms the 
existence of such a transition. 

The adsorption-desorption transition in the model is 
a classical second-order transition. One therefore ex- 
pects a Landau type theory to be applicable for the CEP. 
A CEP is described by a eighth-order Landau function 
F = t(j> 2 + <fi 4 + wcj) 6 + <f> s in terms of a suitable order 
parameter [1]. The first order transition takes place be- 
tween two ordered phases (0 ^ 0) with w > 0, t < or 
between an ordered phase (<fi ^ 0) and a disordered phase 
for w < 0, t > 0. A second order transition takes place be- 
tween the other ordered phase and the disordered phase at 
t = 0. The CEP is at t = 0, w = -3/V2. The first order 
line develops a singularity associated with the behaviour 
of the specific heat across the critical line. Even though 
the first order lines are continuous with the same slope, 
the curvatures are different because the specific heat has 
a jump discontinuity across the critical line. 



By fitting the first order boundaries near the CEP at 
(T e ,£fe) (for e w = 1.8), we find 

gttU (T) = a(T - T e ) + b t , u (T - T c ) 2 + . . . , (7) 

with a = -0.87 ± 0.01, b e = -0.49 ± 0.03, b u = 0.024 ± 
0.019, where the subscript u is for T > T c and I for T < 
T c . This jump in the second-order is consistent with the 
prediction of the Landau theory. Fig. 5 shows the fits 
near the critical end point. 

To summarize, we established the possibility of four dif- 
ferent phases in a set up to unzip an adsorbed dsDNA by 
pulling a single strand. We find that depending upon the 
relative strengths of the binding on the surface e w and the 
pairing energy e\,, either all the four phases or a few of 
them, are present in the phase diagram. By keeping e b 
fixed, we find that a critical end point is present in the 
phase diagram for a wide range of e w . Furthermore, for 
a narrow range of e w , we also have a triple point in the 
phase diagram. It seems that the unzipping of an adsorbed 
dsDNA by pulling a single strand can be a potential can- 
didate to explore the critical end point and this will open 
up a new vista for single molecular spectroscopy. 

Note Added in Proofs: In a recent paper Marenduzzo 
et al. [36] also obtained eq. (5) in a different context of 
melting of a stretched DNA. 
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Fig. 1: Schematic diagram of DNA adsorbed on the surface 
(shaded region). One end of the DNA is always kept anchored 
on the surface at the origin. The free strand (denoted by 1) 
can gam energy Cw for every contact with the surface (i.e. 
x\ — 0). An external force g (shown by arrow) is applied at 
the free end of the pulled strand (denoted by 2 and shifted by a 
unit distance to make it visible). The bold dotted lines denote 
the base pairing (energy — et,) between the two strands of the 
DNA. For all figures, we take ks = 1 and e& = 1. 
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Fig. 2: g vs T phase diagrams (a) for e w = 1 and (b) for 
e w — 1.8. Co represents the critical end point. The reentrance 
on the phase boundary separating phases Z a and Zd is shown 
in the inset. The points are from the transfer matrix and lines 
are Eqs. (4) and (5). 
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Fig. 3: (a) The nature of data collapse of the extensibility at 
T = 1.5 for TV = 1000, 2000 and 3000 at transition Sz. (b) 
T u vs e w curve from data collapse (points) and the exact curve 
(solid line) from Eq. (6). 
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Fig. 4: The {x t )/N vs g (i = 1,2) isotherms for TV = 1000 (a) 
for the 1 + 1 dimensional case at temperatures, T — 0.5 and 
1.5 for e w = 1.8. (b) for the 2 + 1 dimensional case. In both (a) 
and (b), the big squares (for x{) and circles (for X2) show the 
averages obtained by Monte Carlo simulations and the upper 
and the lower triangles are the estimates given by the multiple 
histogram technique at various g. (c) {x\$)/N vs T for the 
DNA of length N = 3000 at g = 0.925. The solid (dashed) line 
represents the free (pulled) strand. The inset shows {xi)/yN 
vs T for various chain lengths. 



Manipulating a single adsorbed DNA for a critical endpoint 
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Fig. 5: Phase diagram in the neighbourhood of C c to show the 
difference in curvature of the two phase boundaries. The lines 
are the fits of Eq. (7) and the points are from the transfer 
matrix calculations. 
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